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The swimming trajectories of self-propelled organisms or synthetic devices in a viscous 
fluid can be altered by hydrodynamic interactions with nearby boundaries. We explore a 
multipole description of swimming bodies and provide a general framework for studying 
the fluid-mediated modifications to swimming trajectories. The validity of the far-field de- 
scription is probed for a selection of model swimmers of varying geometry and propulsive 
activity by comparison with full numerical simulations. The reduced model is then used 
to deliver simple but accurate predictions of hydrodynamically generated wall attraction 
and pitching dynamics, and may help to explain a number of experimental results. 

1. Introduction 

The swimming kinematics and trajectories of many microorganisms are altered by 
the presence of nearby boundaries, be they solid or deformable, and often in perplexing 
fashion. This activity takes place at extremely low Reynolds numbers, a regime in which 
fluid motion is dominated by viscous dissipation. An important factor for swimming 
at such scales is the long range nature of hydrodynamic interactions, either between 
immersed bodies, or between an immersed body and a surface (Lauga & Powers 2009). 
When the swimming dynamics of an organism vary near such boundaries a question arises 
naturally: is the change in behavior biological, fluid mechanical, or perhaps mediated by 
other physical laws? E. coli cells, for instance, have been observed to swim in large circles 
when in the presence of a solid boundary, which has been accounted for in a purely fluid 
mechanical consideration by Lauga et al. (2006). Other organisms have been shown to 
reverse direction at boundaries by inverting the orientation of flagellar rotation, resulting 
in a departure from the boundary which is clearly not a passive hydrodynamic effect 
(Cisncros et al. 2007). In an attempt to help differentiate such observations, we seek 
a general framework for determining the extent to which fluid mechanics can passively 
alter the swimming trajectories of microorganisms near surfaces. 

Surface effects on motility lead to varied and important consequences in a number of 
engineering and biological systems. Van Loosdrecht et al. (1990) note that surfaces are 
the major site of microbial activity in natural environments, and refer to Harvey & Young 
(1980) who showed that almost all detectable bacteria in a marsh estuary were associated 
with particles. Correspondingly, the attraction of certain microorganisms to surfaces has 
a major impact on the development of biofilms, which can begin with the adhesion of 
individual cells to a surface (Van Loosdrecht et al. 1990; O'Toole 2000). Biofilms are 
responsible for numerous microbial infections, and can play an important role in such 
phenomena as biological fouling (Lynch, Lappin-Scott & Costerton 2003; Harshey 2003; 
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Kolter & Greenberg 2006). A recent review on the mathematical modeling of microbial 
biofilms has been presented by Klapper & Dockery (2010). Meanwhile, in a lab setting 
it is common that microorganisms are in near contact with microscope slides or are 
directed through microfluidic channels in which boundaries can play significant roles. 
The migration of bacteria through small-diameter capillary tubes was studied by Berg 
& Turner (1990), and that of infectious bacteria along medically implanted surfaces was 
considered by Harkes, Dankert & Feijen (1992). More recently, Evans & Lauga (2010) 
have shown that the presence of a wall can lead to a change in the waveform expressed 
by actuated flagella, which in turn results in an increase or decrease in its propulsive 
force depending on the type of actuation. 

The study of microbial attraction to surfaces reaches back to the observations of Roth- 
schild (1963), who measured the distribution of bull spermatozoa swimming between 
two glass plates and found the cell distribution to be nonuniform with the cell density 
strongly increasing near the walls. By modeling swimmers as dipolar pushers, it has been 
argued by Berke et al. (2008) that this hydrodynamic consideration alone can account 
for the attraction. Immersed boundary simulations of swimming bodies with undulating 
flagella have also shown a hydrodynamic attraction towards a wall (Fauci & McDonald 
1995). More recently, Smith et al. (2009) have explored numerically the wall effects on 
geometrically accurate swimming human spermatozoa. They have demonstrated that hy- 
drodynamic interactions can trap the body in a stable orbit near a boundary, in some 
cases with counter-intuitive orientation and at finite separation distance from the wall. 
The numerical results of Shum, Gaffney & Smith (2010) and Giacche, Ishikawa & Yam- 
aguchi (2010) show the existence of a stable swimming distance from the boundary in 
swimming E. coli that depends upon the shape of the cell body and the flagellum. Goto 
et al. (2005) have also detected an equilibrium pitching angle for a given wall separation 
distance. Recent experiments showing the upstream swimming of bacteria in a shear 
flow by Hill et al. (2007) suggest that the geometry and orientation of hydrodynamically 
bound swimming organisms can be important. Meanwhile, Drescher et al. (2011) have 
measured experimentally the flow field generated by the swimming of an individual E. 
coli bacterium near a solid surface and have shown that stcric collisions and near-field 
lubrication forces dominate any long-range fluid dynamical effects on these length scales. 

Other recent studies on the dynamics of swimming bodies near walls of a more theoret- 
ical nature includes work by Zargar, Najafi, & Miri (2009), who studied the dynamical 
motion of a three-sphere swimmer near a wall, and Crowdy & Or (2010), who stud- 
ied a simple two-dimensional model of a swimmer using methods of complex analysis 
(see also Crowdy 2011). A different avenue of inquiry has also seen much recent activ- 
ity, the effect of boundaries on swimming suspensions of microorganisms. For instance, 
Hernandez-Ortiz et al. (2009) have studied model swimmers composed of dipolar pushing 
beads, and have shown that the additional length scales introduced by confinement can 
suppress the onset of large scale structures in the suspension. 

Frequently it is the case that the surface of interest does not impose a no-slip condi- 
tion on the fluid velocity, for instance at a free boundary between water and air. Tuval 
et al. (2005) have considered the development of large scale fluid structures driven by a 
competition between oxygen-taxis near the surface of a sessile drop and gravitational ef- 
fects. Di Leonardo et al. (2011) considered the hydrodynamic interactions of a swimming 
bacterium with a stress-free surface, which can be analyzed by placing a mirror image 
of the swimming organism opposite the free-surface. The circular trajectories studied by 
Lauga et al. (2006) were found to be reversed in this setting. 

Outside of the fundamental benefits of scientific inquiry, a more complete understand- 
ing of the hydrodynamic interactions between self-propelled bodies will continue to drive 
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the development of engineering applications as well. Synthetic swimming particles have 
been designed to perform tasks at an exceptionally small length scales, including chemi- 
cally driven bimetallic nano-rods (Paxton et al. 2004; Fournier-Bidoz et al. 2005; Ruckner 
& Kapral 2007), magnetic nanopropellers (Ghosh & Fischer 2009; Pak et al. 2011), and 
undulatory chains of magnetic colloidal particles (Dreyfus et al. 2005) (see also Wang 
2009). Sorting and rectification devices which lean upon asymmetries in microbial inter- 
actions with walls have been explored by Galajda et al. (2007), while Di Leonardo et al. 
(2010) have considered the driven motion of gear-like ratchets in bacterial suspensions. 
Another application of more recent interest is in the production of biofuels, where sus- 
pensions of algae are shuttled through long channels (see Bees & Croze 2010). Exploring 
the hydrodynamic interactions between self-propelled bodies and surfaces not only allows 
us to understand the biological realm with greater sophistication, but may also allow for 
the development of manmade devices of increasing complexity and creativity. 



In the present study, we utilize a multipole representation of self-propelled organisms 
in order to improve our understanding of swimming behaviors near a surface from a 
generalized perspective. The modeling of swimming organisms by Stokcslet dipolc sin- 
gularities has become commonplace, but here we take one systematic step further in 
the far-field expansion of the flow generated by self-propelling bodies. The inclusion of 
higher order singularities will be shown to have important consequences on swimming 
trajectories. Full scale simulations of the Stokes equations are used as a benchmark to 
explore the regions of validity and limitations of the reduced model for two types of 
model swimmers, namely ellipsoidal Janus particles with prescribed tangential surface 
actuation and bacteria-like spheroid-rod swimmers. The far-field approximation leads to 
very good quantitative agreement with the full simulation results in some cases down 
to a tenth of a body length away from the wall. Exploiting the quantitative predictions 
from our singularity approach, the reduced model is further shown to provide good pre- 
dictive power for the initial attraction/repulsion to the wall and the rotation induced by 
the presence of the wall, and even surface scattering in the particular case of spheroidal 
squirmers. 



This paper is organized as follows. In §11, the Stokeslet and higher order singularity 
solutions of the Stokes equations are introduced, and a general axisymmetric swimmer 
is described in terms of a Stokeslet dipole, a source dipole, a Stokeslet quadrupole, and 
a rotlet dipole. The wall effects on the trajectory of a swimming body are described 
through the contribution of each singularity in §111. In §IV we address the question of 
how accurately this multipole representation of swimming trajectories captures the real 
wall effects observed in full numerical solutions of the Stokes equations. The reduced 
(singularity) model is used to provide a simple description of the wall-induced rotation 
for model Janus particles, as well as to describe the complete swimming dynamics of a 
squirming spheroid in §V. In §VI we consider model polar swimmers that are bacteria-like 
in their geometry, and develop an approximate Faxen Law for their study. In §VII we then 
employ the reduced model to study a transition in the wall-induced rotations experienced 
by bacteria-like swimmers for a critical flagcllum length. We finally discuss the accuracy 
and limitations of the reduced model in describing the geometry and dynamics of trapped 
self-propelled bodies near surfaces. 



4 



S. E. Spagnolie and E. Lauga 



2. Singularity representation of motion in a viscous fluid 

2.1. The Stokes equations and singularity solutions 

The length and velocity scales which describe the locomotion of microorganisms are 
extremely small. The fluid flow generated by their activity is dictated almost entirely by 
viscous dissipation, as summarized by Purcell (1977). The Reynolds number describing 
the flow, Re = p UL/p is likewise very small, where p is the fluid density, fi is the dynamic 
viscosity, and L and U are length and velocity scales characteristic of the organism. The 
swimming of E. coli, for example, is characterized by a Reynolds number Re w 10~ 4 (see 
Childress 1981). The fluid behavior is therefore described well by the Stokes equations, 

V • a = -Vp + fiAu = 0, (2.1) 
V-u = 0, (2.2) 

where a = —p I + 2pE is the Newtonian fluid stress tensor, p is the pressure, u is the 
fluid velocity, / is the identity operator, and E = (Vit + Vw T )/2 is the symmetric rate- 
of-strain tensor. The fluid velocity is assumed to decay in the far-field, and the boundary 
condition assumed on the swimming body depends on the specific organism, as will be 
described below. In situations where we include the presence of a plane wall of infinite 
extent at z = 0, we shall also assume a no-slip condition there (u(z = 0) = 0). Now 
classical treatises on zero Reynolds number fluid dynamics have been written by Happcl 
& Brenner (1965) and Kim & Karrila (1991). 

The linearity of the Stokes equations allows for the introduction and exploitation of 
Green's functions. The description of the fluid behavior far from an actively motile or- 
ganism, for instance, can be described accurately using only the first few terms in a 
multipole expansion of fundamental singularities, which will be our approach here. The 
utilization of fundamental singularities allowed for a series of exact solutions to funda- 
mental problems in Stokesian fluid dynamics to be derived by Chwang & Wu (1975). 

A free-space Green's function for the Stokes equations is derived by placing a point 
force in an otherwise quiescent infinite fluid, fS(x a ) (where 8(x ) is the Dirac delta 
function centered at x ). With the point force directed along the unit-vector e (and 
defining / = /e), the solution to the forced system produces the so-called Stokcslct 
singularity, 

u(x) = J-G(x-x ;e), (2.3) 

where 

^.f 1 / \e ■ (x — x )](x — x )\ . .. 

G(x-x ;e) = -\e+ l —± ^ (2.4) 

is the e-directed Stokeslet, and R = \x — x a \. Derivatives of the Stokeslet singularity 
produce other higher-order singularity solutions of the Stokes equations. The first three 
such singularities are the Stokeslet dipole, quadrupole, and octupole, described by 

Gn(x — x ; d,e) = d ■ \7 G(x — x a ; e), (2.5) 
G Q (x - x Q ; c, d, e) = c • V G D (x - x Q ; <2, e), (2.6) 
G (x - x Q ; b, c, d, e) = b ■ VqGq(x - x a ; c, <2, e), (2.7) 

respectively, where the gradient (Vo) acts on the singularity placement x a . The vectors 
b, c and d indicate the directions along which each derivative is taken. As the solutions 
above are regular outside of the singular point there are many possible identities that 
may be observed by rearranging the order in which these derivatives are taken. Tensorial 
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expressions of the singularities above are provided by Pozrikidis (1992) (see also Chwang 
& Wu 1975), and we have included the full vector expressions of these singularities in 
Appendix A. In addition to the solutions above, there are singular potential flow solutions 
to the Stokes equations which are associated with Laplace's equation (p = in 2.1). The 
source, source dipole, source quadrupole, and source octupole singularity solutions are, 
respectively, 

U{x-x ) = ?ZJ*L, (2.8) 

D(x — x Q ;e) = e ■ VoU(x — x ), (2.9) 
Q(x - x a ; d,e) = d ■ VqD{x - x ; e), (2.10) 
0(x — x ; c,d,e) = c • V a Q{x - x a ; d, e). (2-11) 

The source singularities are related to the force singularities through the relation 

D(x-x o ;e) = -^V 2 G{x-x o ;e), (2.12) 

and its derivatives. A notable combination of the above singularities has been named 
alternately the couplet or rotlet, 

R(x - x D : e) = ^ (G D (x - x ; e^,e^) - G D (x - x D : , e^)) 

= ex{x R r°\ (2.i3) 

where the vectors e, e^e- 11 - form an orthonormal basis with 6X6^= e ±J -. A rotlet 
dipole may then be written simply as 

Rn(x — x a ; d,e) = d ■ WoR(x — x a ; e). (2-14) 

Finally, a combination of the above singularities that we will require later is termed the 
Stresslet, which may be written as 

T(x — x a ; d, e) = 2(d ■ e)U(x — x ) + (Gn(x — x a ; d, e) + Gn(x — x a ; e, d)) (2.15) 
= - 6( " o) [d-(x-x )][e-(x-x )}, (2.16) 

representing the fluid force on a plane with normal d corresponding to the e-directed 
Stokeslet velocity field. 

2.2. Far-field description of a swimming body 

In the present study, we will consider only microorganisms which are to a good approxi- 
mation axisymmetric along a direction indicated by the unit-vector e. Such an organism, 
with its centroid at a point x 01 generates fluid motion in the far-field of the form 

u(x) = a G D (e, e) + [3 D(e) + 7 G Q {e, e, e) + t R D (e, e) + O (\x - x a \- 4 ) . (2.17) 

Here we have introduced the shorthand notation Gd(x — x a ; e, e) = Gt>(e,e). The 
coefficient a has units of [Velocity] x [Length] 2 , while (3, 7 and r have units of [Velocity] x 
[Length] 3 . The values of the coefficients a,/3, 7, and t must be determined for each 
microorganism, and depend on the specific body geometry and propulsive mechanism. 

An illustration of the singularity decomposition above is provided in Fig. 1. At leading 
order, a swimming E. coli organism can be modeled as a force dipole (decaying as 1/R 2 ). 
This leading order representation has been used by many authors to consider the effects 
of nearby walls (Berke et al. 2008), and many-swimmer interaction dynamics (sec for 
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Offset force dipole — > 

force quadrupole 
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Figure f . The fluid velocity far from a swimming E. coli is modeled at leading order as that 
of a Stokeslet dipole. At the next order, the flow in the far-field varies due to the length asym- 
metry between the backward-pushing propeller and the forward-pushing cell body (producing 
a Stokeslet quadrupole), the finite size of the cell body (producing a source dipole), and the 
rotation of the flagellum and counter-rotation of the cell body (producing a rotlet dipole). 



instance Hernandez- Ortiz et al. 2005; Saintillan & Shelley 2008; Hohenegger & Shelley 
2010). A swimmer such as the one illustrated in Fig. 1, in which a flagellar propeller 
pushes a load through the fluid, is generally referred to as a pusher, in contrast to such 
organisms as Chlamydomonas which pulls a cell body through the fluid with a pair of 
flagella. At the next order (decaying as 1/R 3 ), the flow in the far-field varies due to the 
length asymmetry between the backward-pushing propeller and the forward-pushing cell 
body (producing a Stokeslet quadrupole) , due to the finite size of the cell body (produc- 
ing a source dipole) , and due to the rotation of the flagellum and counter-rotation of the 
cell body (producing a rotlet dipole). Vector field cross-sections of the Stokeslet dipole, 
source dipole, and Stokeslet quadrupole singularities are shown in Fig. 2. The strengths 
of these singularities have been measured experimentally for the organisms Volvox carteri 
and Chlamydomonas reinhardtii by Drescher et al. (2010), and for E. coli by Drescher 
et al. (2011). The effects of the Stokeslet quadrupole component of spermatozoan swim- 
ming has been suggested by Smith & Blake (2009), and force-quadrupole hydrodynamic 
interactions of E. coli have been studied by Liao et al. (2011). 

While the flow field is set up instantaneously in Stokes flow upon the variation of an 
organism's geometry, the means of propulsion of a particular organism might be unsteady. 
In general the singularity strengths can be time-dependent, varying for example with the 
different phases of a swimming stroke pattern. As an example, the highly time dependent 
flow field generated by the oscillating motions of C. reinhardtii has been examined by 
Guasto et al. (2010). Nevertheless, for a first broad look at the far-field representation 
above we will restrict our attention to constant values of the singularity strengths for 
the remainder of our study. Also, we have assumed in the description given by (2.17) 
that there are no net body forces or torques on the organism, which would require the 
inclusion of Stokeslet and rotlet singularity terms as well (as explored for the organism 
Volvox by Drescher et al. 2009) . While some organisms are not neutrally buoyant and do 
experience a body force or torque due to gravity, many others (including most bacteria) 
live on such a scale that such effects are negligible. In addition, we assume that there 
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Figure 2. Velocity field cross-sections of (a) a Stokeslet (force) dipole, which decays as 1/R 2 ; 
(b) a source dipole, which decays as 1/R 3 ; and (c) a force quadrupole, which decays as 1/-R 3 , 
all in free-space. Arrow intensity correlates with the magnitude of the velocity. The effects of a 
nearby boundary may be intuited by imagining the wall to follow the streamlines. 



is no mass flux through such mechanisms as fluid extrusion, as studied by Spagnolie & 
Lauga (2010), which can present a source singularity in addition to those included in the 
expression above. 

2.3. The surface effect: Faxen's Law 

In a fluid of infinite extent, the fluid velocity in the far-field generated by an active body 
behaves as described in (2.17). When a boundary such as a plane wall is present, however, 
the velocity everywhere is altered due to the additional boundary condition. Borrowing 
an approach which has seen a long history in electrodynamics, the boundary condition 
on the surface can be satisfied by the placement of additional singularities at the image 
point x* a = x a — 2(x a ■ z)z inside the wall (where z is the unit vector normal to the 
surface). 

The image singularities required to cancel the effects of Stokeslet singularities placed 
parallel or perpendicular to a no-slip wall have been presented by Blake & Chwang (1974), 
each requiring a Stokeslet, Stokes doublet, and source dipole, as described in Appendix 
B. The image system for a "tilted" Stokeslet (a Stokeslet directed at an angle relative 
to the wall) is simply a linear combination of the wall-parallel and wall-perpendicular 
image systems. The image systems for higher order singularities, however, are not simply 
linear combinations of wall-parallel and wall-perpendicular image systems. The images 
for each of the axisymmctric singularities in (2.17) will be denoted by an asterisk (see 
Appendix B). For instance, the effect of a Stokeslet dipole along with its image, evaluated 
on the wall surface z = 0, returns Gd{z = 0; e, e) + G* D (z — 0; e, e) = 0. Likewise, we 
denote by u* (x) the fluid velocity generated by the entire collection of image singularities 
needed to cancel, on the no-slip wall, the swimmer-generated velocity description in 
(2.17). Magnaudet et al. (2003) have recently taken a similar approach to studying the 
deformation and migration of a drop moving near a surface, and have provided a valuable 
review of Faxen's technique. 

The flow generated by the image singularities indicates the alteration to the fluid 
motion everywhere due to the presence of the wall. The effects of this induced fluid 
motion on the swimming trajectory are provided by Faxen's Law, which can be written 
exactly for a prolate ellipsoidal body geometry (Kim & Karrila 1991). For an ellipsoid 
of major and minor axis lengths 2a and 26, with the major axis aligned with the unit 
vector e, the translational velocity u and rotational velocity fl induced on the swimmer 
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due to its experience of the flow u*(x) may be written as 

u = u*(x )+0(a 2 \7 2 u*\ Xo ) , (2.18) 

fi=^Vx u*(x ) + re x (£*(x ) ■ e) + O (a 2 V 2 (V x u*) \ Xo ) , (2.19) 

where x a is the position of the body centroid, E*(x) = (Vm* + [Vtt*] T )/2, and E = 
(1 — e 2 )/(l + e 2 ), with e = 6/a the body aspect-ratio. Denoting the swimming speed 
attained by the organism in free-space by U, we have therefore that the body swims with 
velocity Ue + u and changes swimming direction via e = fl x e. 

At the order of our consideration in this paper (via 2.17) the strength of the singular- 
ities representing the body motion are not changed by the presence of the wall. If the 
singularities required to represent the motion differed in rate of decay by more than one 
degree of separation (we currently include only terms decaying at order R~ 2 and i?~ 3 ), 
then Faxen's Law above would indicate a problematic interaction of the wall effect with 
the measurement of the singularity strengths. The approach above, then, must be han- 
dled with more care in the event that a Stokeslet singularity is required, or if higher order 
terms than those considered here are to be included in (2.17). So long as the distance of 
the body to the wall is sufficiently large relative to the body size, the higher order terms 
in (2. 18)- (2. 19) may be neglected. 

The expressions above do not extend easily to geometries that are not ellipsoidal. In 
order to study a body geometry more like that illustrated in Fig. 1 we will need to 
develop an approximate "Faxen Law." First, however, let us consider the consequences 
of singularity images on a prolate ellipsoidal body to develop some intuition. 



3. The surface effect, singularity by singularity 

A swimming body first begins to experience the hydrodynamic effects of a wall through 
the singularities which decay least rapidly in space. We now list each singularity in 
the multipole representation and describe the corresponding effect on the swimming 
speed and orientation of the body. First, however, the system is made dimensionless by 
scaling velocities upon the free-space swimming speed, U, lengths upon the body semi- 
major axis length a, and forces upon jiaU . Henceforth all variables are understood to 
be dimensionless. The unit vector z is normal to the wall and points into the fluid, and 
the dimensionless distance between the wall and the body centroid is denoted by h (see 
Fig. 3). The pitch angle with respect to the wall is denoted by 6; the body is swimming 
directly away from the wall when 9 — ir/2, directly towards the wall when 9 = —n/2, 
and parallel to the wall when 9 = 0. 

3.1. Force dipole 

Since we will not consider external body forces or torques, the least rapidly decaying sin- 
gularity (in space) generated by the activity of an organism is a Stokeslet dipole directed 
along the swimming direction, Gt>(e, e), which induces a dimensionless attraction to the 
wall (or repulsion from the wall) by (2.18) of the form 

z-u = u z = -^(l-3sm 2 9). (3.1) 

Further details are provided in Appendix B, along with the effects induced not by a 
wall but instead by a stress-free surface such as a fluid-air interface. From (3.1), the 
surface-induced velocity of a pusher (a > 0) is towards the wall when | sin(6 l )| < l/VS- 
For small orientation angles 9^0 (when the body swims almost parallel to the wall, 
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Figure 3. Schematic representation of the generic problem studied in this paper: an inclined 
swimmer near a solid surface. The distance of the body centroid from the wall is denoted by 
h, measured along the direction normal to the wall z. The pitch angle of the body's director 
e with respect to the wall is denoted by 6; the body is swimming directly away from the wall 
when 6 = ir/2, directly towards the wall when 6 = — it/2, and parallel to the wall when 8 = 0. 



and sin(6*) ~ 8), combining the wall effect above with the vertical component of the free- 
space swimming speed, sin(#), we see that the body will move towards the wall when 
8 < 3a/ (8h 2 ). Hence, for a pusher that is swimming nearly parallel to the wall, the first 
effect of the hydrodynamic interaction with the wall is an attraction. It has been argued 
by Berke et al. (2008) that this hydrodynamic consideration can account for observations 
of the entrapment of E. coli near surfaces, as well as the observations of Rothschild 
(1963), who measured the distribution of bull spermatozoa swimming between two glass 
plates and found the cell density to increase near the walls. Experimental measurements 
of force dipole strengths generated by swimming E. coli were evaluated by Drescher 
et al. (2011), who found rotational diffusion to dominate hydrodynamic effects in that 
particular regime. 

Now, what of the body orientation dynamics? The pitch angle 8 is assumed to be 
constant in time absent the presence of a wall, so the variation in 8 is due only to wall- 
induced rotational effects. The leading order effect is again that of the slowly decaying 
Stokeslet dipole term, which generates the rotation rate 



where the approximation given is appropriate when the body is nearly parallel to the 
wall (see Appendix B for the full expression). Hence, for 8 ~ the induced rotation 
acts to align the body with the wall for a > (pushers), and perpendicular to the wall 
for a < (pullers), with no qualitative dependence upon the aspect-ratio of the body. 
The rotation induced by the force dipole may be intuited based on a consideration of 
the velocity fields shown in Fig. 2a; imagining a wall to follow the streamlines, the body 
is seen to be drawn into the wall, and based on velocity gradients to rotate towards 
8 = 0. The nature of the wall effect is even more predictable for the source dipole and 
force quadrupole from a similar consideration of Figs. 2b-c, as we will show. As we now 
proceed to consider the next order of singularities, we shall find that the leading order 
wall effects described above can be rather deceptive if they are used to predict the full 
trajectory of a given swimmer. 




(3.2) 
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3.2. Source dipole 

As a swimming body comes into closer contact with a wall, or when the body is swimming 
parallel to the wall, higher order singularities will begin to affect the trajectory of the 
swimming organism. The source dipole singularity, which enters due to the presence of 
a cell body of finite size (such as the spherical head in Fig. 1) contributes an induced 
attraction/repulsion relative to the wall of the form 

For a swimmer such as that shown in Fig. 1 it is common to have /3 < 0, since motion 
of an inert sphere through a viscous fluid can be represented by a Stokeslet singularity 
and a source dipole with j3 < placed at its center (Kim & Karrila 1991). In this case, 
the source dipole term contributes a wall repulsion when the body is pitched away from 
the wall ( "nose up" ) , and contributes an attraction to the wall when the body is pitched 
down ("nose down"). 

When 8 = 0, the Stokeslet dipole term no longer determines the rate of rotation; 
instead, it is set by the higher order singularity structure. The rotational velocity induced 
(for 9 ~ 0) by the source dipole is 

3/3 / 3r 



9 ~w\ l + ^)> (3 - 4) 

which acts to rotate the nose downward towards the wall if j3 < 0. In general, these 
effects will compete with those generated by the Stokeslet dipole and quadrupole to set 
the trajectory and equilibrium states of the self-propelled body. 

If instead a body swims by activity on its surface (as is the case for ciliated organisms, 
often modeled as so-called squirmers), we can find [3 > 0, as will be shown in the following 
section. The noted effects above are thus reversed for such an organism, which we will 
explore in greater detail in §V. 

3.3. Force quadrupole 

At the same order of decay as the source dipole, the Stokeslet quadrupole enters and 
induces a wall-perpendicular velocity of the form 

^=7^(7-9sin 2 fl), (3-5) 
and contributes an induced rotation rate of 



for 8^0. The attraction/repulsion and induced rotation rate depend on the sign of 
7, which itself depends upon the fore-aft body asymmetry, indicated in Fig. 1. From 
studying swimmers with exact singularity expressions (to be described below), we expect 
7 < for such swimmers as shown in Fig. 1 with large cell bodies and short flagella, 
and 7 > for those with small cell bodies and long flagella. Like the source dipole, this 
singularity also acts to rotate the swimmer when 8 = 0. 

3.4. Rotlet dipole 

The rotlet dipole term can account for at least one surprising behavior of locomotion 
near surfaces, the circular swimming trajectories of E. coli as studied by Lauga et al. 
(2006). For a body swimming parallel to the wall, 8 = 0, the rotation about the z axis 
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is given by 

the effect disappearing for infinitely slender swimmers (for fixed r). Fixing the distance 
to the wall h, the body is thus predicted to swim in circles with a (dimensionless) radius 

* " sFRTTo- (3 ' 8) 

The cell bodies of E. coli bacteria rotate clockwise as seen from the distal end (behind 
the organisms) during their forward swimming runs, and the net torque is balanced by 
the counter-clockwise rotation of the propelling flagellar bundle. This situation corre- 
sponds to r > 0, and hence (3.7) predicts a large clockwise circular trajectory (as seen 
from above) parallel to the plane of the wall, which is consistent with the experimental 
observations of Lauga et al. (2006) . 

Note that the same organism moving near a stress-free boundary (like an air-water 
interface) experiences a passive rotation in the opposite direction (see Appendix B), 
as studied by Di Leonardo et al. (2011). While the rotlct dipolc contributes to three 
dimensional swimming dynamics, this component of the propulsion has no bearing on 
the wall- attraction/repulsion or pitching dynamics of a swimmer: u = and 9 = 0. 
This is assured by the kinematic reversibility of Stokes flow. We will focus here on wall- 
attraction/repulsion and pitching dynamics, and thus for the remainder of our consider- 
ation we will set r = in (2.17). 



4. Where is the multipole singularity representation valid? 

The central question that we wish to answer in this paper is: how accurately are the wall 
effects predicted by the multipole singularity representation (2.17) and accompanying 
Faxen Law (2.18)-(2.19)? In order to provide an answer, we must have at our disposal 
a means of computing the full fluid-body interaction. For this purpose we will utilize 
the method of images with regularized Stokeslets, as derived by Ainley et al. (2008), 
in which the boundary integral formulation of the Stokes equations is accompanied by 
image singularity kernels which cancel the fluid velocity on the wall. In this approach the 
wall need not be discretized. We will make a necessary adjustment to this framework to 
allow for the inclusion of a slip velocity on the body surface. Equipped with a means of 
computing solutions to the full Stokes equations in the half-space bounded by a no-slip 
wall, we can study the accuracy of the multipole representation and wall effects for a 
selection of model swimmers. 

4.1. Model swimmers 

In order to probe the accuracy of the far-field representation we will enlist the help 
of model "Janus swimmers," as illustrated in Fig. 4. The model swimmers are prolate 
ellipsoids, chosen so that the Faxen expressions (2.18)-(2.19) may be applied exactly. For a 
given pitch angle 9, the surface is parameterized as x = cc + cos(£)e + esin(£;)[cos(ri)e- L + 
sin(?7)e- L - L ], with £ e [0, 7r] and n G [0, 2tt], and where e, e , c 11 - form an orthonormal 
basis with 6X6^= e ±J -. Recall that e is the body aspect-ratio. The unit tangent vector 
everywhere on the surface is denoted by s(£, 77), and the unit normal is denoted by n(£, 77). 

The propulsive mechanism is a prescribed axisymmetric distribution of a slip velocity 
u s which acts tangentially to the body surface for — 1 < cos(£) < £, where I £ [—1, 1] is a 
dimensionless "activity length," while the remainder of the body surface is inert (where 
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2e 




£ = 




Figure 4. A selection of model swimmers are illustrated. Spheres (e = 1) and ellipsoids (e = 1/4) 
are shown with activity lengths £ = —1/2,0, 1/2, and 1. Each body shown shuttles fluid along 
the active part of the surface to the left (the propulsive activity is indicated by arrows), and 
thus swims to the right. The gray regions are inert, where a no-slip condition is assumed. 



a no-slip condition is applied). The entire body surface is active for £ = 1 (a squirmer). 
More specifically, the prescribed slip velocity distribution is chosen to be 



1 2(l+cos£)(£-cos0 
(l + e 2 )-(l-e 2 )cos(20 



(4.1) 



The constant u® is selected so that the free-space (no wall) swimming velocity is unity, 
and is determined numerically. A no-penetration condition is applied on the entire body 
surface. 

The squirmer model of ciliated organisms was introduced by Lighthill (1952) and ex- 
tended by Blake (1971). Squirmer models (where either the slip velocity or the surface 
stress is specified) have been used to study multiple-organism interactions by Ishikawa, 
Simmonds, & Pedley (2006) and Kanevsky, Shelley, & Tornberg (2010), hydrodynam- 
ically bound states by Drescher et al. (2009), efficiency optimization in ciliary beating 
by Michelin & Lauga (2010), fluid stirring effects by Lin, Thiffeault, & Childress (2011), 
and motion in a polymeric fluid by Zhu et al. (2011). Swimmers with partially activated 
surfaces {£ < 1) have recently been designed and studied with great enthusiasm; see for 
instance the work of Paxton et al. (2004), Golestanian, Liverpool, & Ajdari (2007) and 
Jiang, Yoshinaga, & Sano (2010), where the activity is generated by self-phoretic and 
thermophoretic surface effects. Migration of similar 'slip-stick' spheres in an ambient flow 
has been studied by Swan & Khair (2008). 



4.2. Full numerical simulation approach 

An application of Green's theorem to the Stokes equations (2.1) reveals a representation 
of the fluid velocity everywhere based solely on integrations of the stress and velocity on 
the immersed boundaries (the swimming body, in this case) (Pozrikidis 1992). Accounting 
for the presence of the wall by including image singularities, the fluid velocity everywhere 
may be written as 



u(x) = JC[f](x)+)C*lf](x) - {T[u]{x) + T* [«](*)) , (4.2) 
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where 



£[/](*) = ^ / G{x-y;f{y))dS v , (4.3) 



IdD 

1 



8tt 



£*[/](*) = —/ G*(x-y;f(y))dS y , (4.4) 



T[u]{x) = ^-{ T(x-y;n(y),u{y)-u(x))dS v , (4.5) 

T*[u](x) = ^-[ T*(x-y;n(y),u(y)-u(x))dS y , (4.6) 

with / the dimensionless fluid force per unit area, G the Stokeslet singularity (2.4), 
T the Stresslet singularity (2.15), and dS y the differential surface area element for the 
integration variable y. The image singularities Gd* and T* ensure that the no-slip con- 
dition is satisfied on the wall at z = 0, and are provided in Appendix B. The single-layer 
integral JC is weakly singular, which presents both theoretical and numerical difficulties. 
One approach to computing this integral is through the use of a regularized kernel, G$, 
where a small regularization parameter 5 is introduced. This is the approach of Ainlcy 
et al. (2008), who derive the necessary adjustments which must be made to account for 
this regularization, and who also derive the image-singularities (JC*) which must accom- 
pany such an approach when a wall is present. This is the approach taken in the present 
work, though we must include the double layer integral T and its image in order to 
accommodate the slip velocity u s . 

We briefly recount the method of images with regularized Stokeslets as presented by 
Ainlcy et al. (2008) (itself a modification to the method derived by Cortez 2002). The 
surface is discretized by M points, located at x k o for k — 1,2, M. For a given point x in 
the fluid or on the body or wall, we define x* k = x — Xkfi, and Xk = x — (x k fi — 2(z-x k fi)z). 
Absorbing the surface integration into the force so that we may simply write j gD f dS y = 
X)fc /&! an d choosing a blob function of <f>(x) = 15e 4 /87r[r 2 + e 2 ] 7 / 2 (used to spread the 
singular effect of a point force to a small finite area) , it may be shown that 

M 

£[/](*) = £ \hH x (\xl\) + (f k ■ x* k )x* k H 2 (\x* k \)} , (4.7) 

k=l 

£*[f]( x ) = - ifkHi(\x k \) + (f k ■ x k )x k H 2 (\x k \)} - /ifc[t/ fe £>i(|x fe |) + (g k ■ x k )x k D 2 {\x k \) 



-2h k 



+ H 2 (\x k \ 



(L k x x k ) + 2h k (g k ■ z)x k H 2 (\x k \) + (x k ■ z)H 2 (\x k \ 



x H'A\x k \) , ^, x H'(\x k \)i , 

+ (g k ■ x k )z ) Kl , 17 + (x k ■ z)(g k ■ x k )x k ™ * 1; , 4.8 
\x k \ \x k \ 1 

where S is a regularization parameter (discussed below), g k — 2(f k -z)z — f k , L k = f k xz, 
and 

1 5 2 1 

ffl(r) = fa^FW 2 + Sf^FW 2 ' = 8^(r 2 + 5 2 )3/ 2 ' (4 ' 9) 

1 3(5^ 3 

Dl W = 4^/i(r 2 + <5 2 ) 3 / 2 ~ 4^/i(r 2 + <5 2 ) 5 /2 ' = ~ 4 ^( r 2 + ,52)5/2 ' (4 - l ()) 

Having subtracted off the velocity at the target point x, the integrand of the double 
layer integral T[u] is finite with a jump discontinuity at y = x (Power & Miranda 1987; 
Pozrikidis 1992). Insertion of rigid body motion velocities into the integrals returns zero, 
so only the tangential slip velocity u s need be considered. Since the integrands in T[u] 
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and T* [u] are known and finite, they are computed using adaptive quadrature to com- 
puter precision accuracy. 

Having absorbed the surface integration into the definition of fk, the net (dimension- 
less) force on a boundary is computed simply as 



M 

f( X )ds = y2fk, (4.H) 

dD 



fe=l 

which must return zero in the case of self-propelled swimming (where body forces such 
as gravity have been neglected, as previously noted). 

For a given body position and orientation, a linear system must be solved to determine 
the swimming velocity and rotation rate, U and fl, along with the scaled force fk, via 
the boundary integral relation in (4.2). The linear system is closed by requiring that the 
boundary conditions hold as follows. Denoting the inert part of the body by dDi and 
the active part of the body by 3D a, 

u(x e <9L>/) = U+ fl x (x -x a ), (4.12) 
u(x e dD A ) = U + fl x (x - x a ) + u B (£, rj). (4.13) 

The no-slip condition on the wall is satisfied automatically by the inclusion of the image 
kernels. Continuing to follow Ainley et al. (2008), the ellipsoidal surface is discretized 
by dividing the azimuthal angle £ into N% + 1 points, £ TO = irm/N^ for m = 0, 1, N%. 
At each station £ TO , the polar angle is discretized into N m points, i] n = 2im/N m for 
n = 1,2, ...,N m , where N m is the smallest integer larger than 2A / £sin(£ m ). Taking as a 
representative discretization size hd — \J S/N', with N' the total number of gridpoints 
and 5* the (dimcnsionless) ellipsoidal surface area, the regularization parameter chosen 
for all the problems considered herein is 6 = 0.8h d 9 , and we set N% = 30 to capture the 
swimming behavior with sufficient accuracy for our purpose. 

4.3. Computing the singularity strengths: a, (3, and 7 

In order to compare the full system with the singularity representation in (2.17) we 
must determine the singularity strengths for the swimming bodies under consideration. 
Having recovered f(x) numerically, these singularity strengths may be computed as 
follows. Assuming that the body is directed along e — x and centered about x = 0, for 
the sake of presentation, we have by an involved inspection (expanding (4.8) for \x\ » |0| 
in free-space, and matching the x 3 term for the dipole and the x 4 term and (y 2 + z 2 )x 2 
terms for the higher order singularities): 

a = [ (3xf x -x-f)dS+^-f n x u x dS, (4.14) 

lO 71 " JdD 87r JdD 

„2 l™|2\ 



P = j7T {x z -\x\')f x dS 

ID 71 " JdD 



-l 

6ti" Jdl 

— I 2x n x u x — {u x x ■ n + n x x ■ u — 2x n ■ u) dS, (4-15) 

i7r JdD 



1 

327 



(5x' -\x\ z )f x -2x (x-f)dS 

dD 



■—— { hxn x u x — [u x x ■ n + n x x ■ u + xn ■ u) dS. 

871" JdD 



(4.16) 



For purely rigid body motion, u = U + fl x (x — x a ), the integrals involving u vanish, 
so we need only insert the slip velocity u s into these integral expressions. Note that the 
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Figure 5. The singularity strengths a,/3, and 7 are shown for the model swimmers illustrated 
in Fig. 4, but for a wide range of aspect-ratios e and activity lengths I. The Stokeslet dipole 
strength, a, is positive for all such swimmers (they are all pushers). The source dipole strength, 
j3, and the Stokeslet quadrupole strength, 7, change sign where indicated by dashed lines. 



expressions above are geometry dependent; other expressions would need to be derived 
for non-ellipsoidal body shapes. 

The singularity strengths ct,/3, and 7 computed for a range of aspect-ratios e and 
activity lengths £ are shown in Fig. 5. For all such model swimmers as shown in Fig. 4, 
we find that a > (the swimmers are pushers). However, both j3 and 7 change sign, 
as indicated by dashed lines. For a body that is primarily inactive (i w — 1), the source 
dipole term is negative, (3 < 0, which can be predicted by considering the exact singularity 
solution for a sedimenting solid ellipsoid as derived by Chwang & Wu (1975). However, 
when the body is primarily active (£ « 1), we find /3 > 0, which can be predicted from 
the exact solution for a squirming ellipsoid as derived by Keller & Wu (1975). The source 
dipole term decreases in magnitude with e 2 , and also decreases in magnitude as the body 
becomes less geometrically active (£ — > —1). Meanwhile, 7 is negative for bodies with 
small active surface areas £ < and positive for bodies with large active surfaces £ > 0. 

As a simple example, consider a slender rod of dimensionless length 2 which satisfies 
a no-slip condition for arc-lengths s € [£, 1] and has specified active forcing e • / = 
— J-"/(l + £) for arc-lengths s e [— 1,1]. A leading order approximation (in the small 
aspect-ratio of the body) of the force on the no-slip part of the body is simply F /{!—£). 
From (4.14)-(4.16) we find that a = J7(8tt), f3 = 0, and 7 = £T/(24ir). The swimmer 
is a pusher for T > 0, has no source dipole term due to its slenderness, and generates 
a Stokeslet quadrupole term which changes sign precisely when £ = 0, distinguishing 
bodies which are more active than inert or vice versa. The direction of rotation induced 
on such a swimmer by the presence of the wall thus depends here critically upon the 
asymmetry of the propulsive mechanism [through (3.6)], and not the body geometry, 
which is symmetric about its centroid. Namely, when the swimmer is parallel to the wall 
(0 = 0) then 9 > for T£ < and 9 < for T£ > 0. A long inert segment pushed by 
a short active segment (J 7 > 0) leads to an upward pitching motion, while a short inert 
segment pushed by a long active segment leads to a downward pitching motion, with the 
situation reversed if the swimmer is a puller ( T < 0) . 

The singularity strengths for a squirming (potential flow, £ = 1) spheroid are even 
simpler, with a = 7 = 0, and f3 > 0. Keller & Wu (1975) have shown that such a 
swimmer may be described in free-space exactly as an integration over source doublet 
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singularities distributed along the body centerline. Specifically, writing 

u(x) = ^J {c 2 - S 2 )D(x - S e)ds = /3D(x ) + 0(h- 5 ), (4.17) 

2e 2 c 3 

13 = 6c-3e 2 log([l + c]/[l -c])' (4 ' 18) 

where c = \/I — e 2 is the ellipsoidal focal length, then on the body boundary the velocity 
satisfies the rigid body motion (with speed unity) and the tangential slip velocity pre- 
scribed in (4.1) with I = 1. For a spherical squirmer we have = 1/2, and u° = 3/2 in 
(4.1), and the 0(h~ 5 ) term above is precisely zero. A squirming potential flow spheroid 
therefore always rotates away from the wall surface [through (3.4)], which we will explore 
in greater detail in §V. 



4.4. How accurate is the far-field representation? 

Given the computed singularity strengths for the model Janus swimmers we can compare 
the far-field representation of the swimming dynamics with the full simulation results. 
Figure 6 shows the "horizontal" (wall-parallel, in the direction x) and "vertical" (wall- 
perpendicular, in the direction z) velocities and the rotational velocities of the eight 
model swimmers shown in Fig. 4, for a range of distances from the wall, and with pitching 
angle fixed at 9 = it/ 6. Here h w is the minimal distance of the centroid to the wall without 
penetration, h w = \/(l + e 2 ) — (1 — e 2 ) cos(20)/\/2- The far-field approximation from 
(2.17) generally matches the full simulation results with great accuracy in this range. 
Inspection of log- log plots (not shown here) affirm that the singularity strengths have been 
computed correctly, or alternatively verify the accuracy of the full numerical simulations. 
All the bodies with the exception of the squirmers (for which a = 0) swim horizontally 
with greater velocity due to the presence of the wall at this pitch angle. Also note that 
the spherical squirmer is the only swimmer of the eight that rotates away from the wall 
at this pitching angle; the rest all passively rotate towards the wall. 

Importantly, we find that the rotation rates are very small relative to the swimming 
speed, even when the swimmer is as close to the wall as one body length away. For 
swimming organisms directed towards the wall a collision may be inevitable. Meanwhile, 
when a body is swimming parallel to a wall, this small rotation and trajectory adjustment 
may be sufficient to entirely determine whether or not the swimmer collides with the wall 
and enters a hydrodynamically or otherwise bound state, or possibly swims away. 

How close to the wall can the far-field expressions be extended? Figure 7 shows the same 
measurements but for a range of even smaller body distances from the wall. We find that 
the far-field approximation is surprisingly accurate, in some cases all the way down to a 
distance of one-tenth a body length from the wall. In particular, the horizontal swimming 
velocities of both the spherical and ellipsoidal swimmers is generally matched quite well, 
as are the vertical and rotational velocities of the spherical bodies. The behaviors of the 
ellipsoidal bodies are still captured qualitatively in this regime, but we begin to see a 
departure of the far-field approximation from the full simulation results. This is to be 
expected, as the far-field approximation does not capture effects such as lubrication, and 
higher order terms in the multipole expansion eventually become comparable to those 
kept in (2.17). For the swimmers considered here, however, the far-field approximations 
appear to be surprisingly accurate, and we expect that much can be predicted using 
the simple framework encapsulated by (2.17). This said, the predictions can become 
qualitatively inaccurate when the body is almost touching the wall, yielding non-physical 
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Figure 6. The wall-parallel ("horizontal," U x ) and wall-perpendicular ("vertical," U z ) velocities 
and rotational velocities (8) of the eight model swimmers shown in Fig. 4 are shown for a range 
of distances to the wall h, with computed values indicated by symbols and far-field predictions 
indicated by lines. The pitching angle is fixed to 9 = n/6. Top row: spherical swimmers. Bottom 
row: spheroidal swimmers. The far- field approximation is seen to reproduce the full simulation 
results, and generally remains valid for the range considered. The scaling h w is the smallest 
value of h for which there is no wall contact (defined in the text). 



results in some cases, so care clearly must be taken when attempting to apply the far-field 
theory to near-wall or near-swimmer interactions. 



5. Employment of the reduced model: Janus swimmers 

5.1. Equilibrium pitching angles 

Our interest now turns to the rotations induced on a self-propelled body by the presence 
of a wall. Absent any boundaries, the swimmer is assumed to move along a straight 
path, with 8 = 0. When a wall is present, the body will pitch up away from the wall 
or pitch down towards the wall due to hydrodynamic interactions. For a given fixed 
distance from the wall, h, this rotation will persist until an equilibrium angle 6 is reached. 
Stable equilibrium angles so found using the full numerical simulations are shown in 
Fig. 8a, where we have fixed the centroid distance to one body length away from the 
wall, h = 2. Figure 8b shows the same plot but generated using the predictions of the 
far-field theory. Both the full simulations and reduced model predict that there exists a 
class of swimming bodies for which the hydrodynamic interaction with the wall results in 
a negative equilibrium pitching angle. Namely, those that are simultaneously slender and 
more geometrically active than inert, though only to a point. For potential flow squirmers, 
£ — 1, we find that for any aspect-ratio the pitching equilibrium angle is oriented directly 
away from the wall. Inclusive of the potential flow squirmers, we detect a barrier above 
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Figure 7. The horizontal, vertical, and rotational velocities of the eight model swimmers shown 
in Fig. 4, just as in Fig. 6, but for smaller distances from the wall. Computed values are indicated 
by symbols and far-field predictions by lines. The far-field approximation is surprisingly accurate, 
in some cases all the way down to one-tenth the body radius distance from the wall. Other 
velocities begin to vary not only quantitatively but qualitatively as the body approaches wall 
contact. 



which the only equilibrium angle is 9 = 7t/2, and we note that the transition from a 
small orientation angle to no equilibrium other than 9 — n/2 is extremely sensitive to the 
geometry of activation along this curve (a small change in I can produce a dramatically 
different behavior). 

In order to derive a simple analytical estimate of the equilibrium angle, we linearize 
the rotational contribution of each singularity about 9 = 0, yielding 

and hence an equilibrium swimming angle 9 (for a fixed distance h) is predicted to exist 
at 

^ akdlr/2) (7-l+j(U7-6fl). (5.2) 

The estimate from the linearized far-field theory is show in Fig. 8c. The results match 
the full far-field theory and simulations well where the equilibrium angle is small, as 
expected, showing an overestimate in the angle for bodies which are less geometrically 
active (£ ~ —1). That the completely active bodies (£ = 1) have no stable equilibrium 
pitching angle outside of 9 = 7r/2, in which the body is swimming directly away from the 
wall, is explored in greater detail in the following section. 

The dashed curve marking the transition from pitched-up to pitched-down equilibria 
also separates bodies that will rotate away from the wall and those that will rotate 
towards the wall when 9 = (which may be inferred by continuity) . Pausing to consider 
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Figure 8. Comparison of the full simulations to the far- field predictions for fixed distance 
h = 2. (a) Computed contours of the stable equilibrium pitching angle at fixed distance h — 2. 
(b) The same, as predicted with the far-field theory, (c) Analytically predicted equilibrium angle 
using the far-field approximation, linearized about 8 = (taking the minimum value between 
the prediction in (5.2) and 7r/2). Bodies that are both slender and sufficiently active (but not 
completely active, I = 1) exhibit pitching equilibria with their noses down towards the wall. 
Bodies that are not sufficiently slender or not sufficiently active exhibit pitching equilibria with 
their noses turned up away from the wall. 



this rotation rate for wall-parallel swimming (6 = 0), we refer to Fig. 5 and note that 
the rotation is towards the wall in a region where < j3 -C 1 and 7 > 0, and also where 
r = (1 — e 2 )/(l + e 2 ) > 0.5. For i > 0, Fig. 5 shows that j3 grows more rapidly than 7 
as the aspect-ratio is increased. For £ > 0, then, both (7 — f3) and r(llj — 6/3) /i will 
become negative for sufficiently large aspect-ratio, and 9 will become positive [see (5.1)]. 
In other words, the effect of the source dipole will overwhelm the effects of the Stokeslet 
quadrupole for sufficiently large aspect-ratio when I > 0. When I < the situation is 
reversed: 7 decreases more rapidly than does (3 as the aspect-ratio is increased. In this 
regime the Stokeslet quadrupole effects dominate those of the source dipole, both (7 — 0) 
and -T(ll7 — 6/3)/4 are negative in (5.1), and thus 6 > 0. 

To summarize, bodies that are both slender and sufficiently active exhibit pitching 
equilibria with their noses down towards the wall, though completely active bodies [i = 1) 
always rotate away from the wall. Bodies that are not sufficiently slender or are not 
sufficiently active exhibit pitching equilibria with their noses turned up away from the 
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surface. In addition, there is a boundary in parameter space beyond which the only stable 
equilibrium angle is 9 = tt/2, with the body swimming directly away from the wall. 

5.2. Full swimming trajectories of squirmers 

Unsurprisingly, the far-field approximation cannot generally be counted on for quantita- 
tive (and in many cases qualitative) predictions of the entire swimming behavior when 
the body is in very near contact with the boundary. The exact form of the propulsive 
mechanism (in this case the form of the tangential slip velocity) will specify the nature 
of near-wall contact, be that a hydrodynamically trapped state or a trajectory that leads 
to a departure from the surface. However, we have found one class of swimmers for which 
the far-field approximation can be used to predict the entire interaction with the bound- 
ary, namely for squirmers (£ = 1), which we now describe in detail. That the far-field 
theory provides an accurate depiction of the full dynamics of a treadmilling swimmer 
was found in a two dimensional setting by Crowdy & Or (2010) and Crowdy (2011). 
The interaction of a squirmer with a wall has also been studied recently by Llopis & 
Pagonabarraga (2010). 

We comment briefly on the numerical method. Time does not enter into the Stokes 
equations explicitly, and since the means of propulsion studied here is steady there are 
no variations in the dynamics with time outside of the trajectories described by the 
distance of the centroid from the wall, h(t), and the pitching angle, 9{t). An adaptive 
time-stepping algorithm for stiff systems (odel5s in Matlab) is used to integrate the 
swimming trajectory to small enough error tolerance so that the sole error in the dynamics 
is due to discretization errors and the associated quadrature errors in evaluating JC[f] 
(through the regularization parameter S). Hydrodynamic interactions with the wall are 
sufficient to prevent body-wall collisions in some but not in all cases. Following Brady 
& Bossis (1985) and Ishikawa & Pedley (2007), we include a screened electrostatic-type 
body repulsion force which acts only at very small distances from the boundary of the 
form F re .p = Ae Bd /[l - e Bd ]z, where d is the minimum distance between the body 
surface and the wall, and we take A = 100, B = 10. These values are selected so that the 
body does not come closer than approximately h/h w = 1.05, where the numerical method 
for computing the fluid velocity just begins to lose accuracy (though not dramatically; 
see Ainley et al. 2008). For elongated squirmers (e < 1) an associated torque is included 
as well. More physically realistic near-contact interaction effects have been discussed by 
Poortinga et al. (2002). 

We first describe the results of the full simulations by focusing on a spherical squirmer 
(e = 1). By scanning the parameter space of distances h and pitching angles 9, we have 
observed that the body rotates away from the wall regardless of its distance from the 
boundary and orientation (save for the special case of swimming directly towards the wall 
at 9 = —tt/2, though this orientation is found to be unstable). Setting the body initially 
at a distance one body radius from the wall (h = 2), for initial pitching angles 9 > —0.4 
the body rotates as it moves through the fluid and docs not come into close contact 
with the surface. It then settles into a final pitching angle 9f > as it swims away from 
the wall, never to return. This scattering angle, as determined from full simulations, is 
shown in Fig. 9a as a solid line. For initial angles 9 < the body swims towards the wall, 
which increases the rotational effect on the body and hence increases the final pitching 
angle once it departs, accounting for the non-monotonicity of the scattering curve. Three 
trajectories have been picked out with the intention of illustrating this non-monotonicity 
and are included as Fig. 9b. 

Meanwhile, for initial pitching angles 9o ^ —0.4 the squirmer "impacts" the wall (real- 
ized here by blocks of time in which there is negligible wall-perpendicular velocity), but 
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Figure 9. (a) The scattering angle exhibited by a squirming sphere (I = 1 and e = 1), for an 
initial centroid distance h = 2. Solid line: results from the full simulation. Circles: analytical 
far-field prediction, neglecting the surface attraction/repulsion, from (5.7) and (5.10). (b) The 
trajectories of spherical squirmers with three different initial pitching angles, illustrating the 
non-monotonicity indicated in (a). A trajectory that takes the body nearer to the wall leads to 
a greater net rotation before the swimmer escapes. 

continues to rotate while in near-wall contact. Eventually the squirmer escapes the sur- 
face and swims away until settling to a final pitching angle Of s=s 0.4. For non-impacting 
trajectories, the swimming trajectory must be symmetric about the point of parallel 
swimming, = 0, by the time-reversibility of the Stokes equations (reversing the direc- 
tion of time is indistinguishable from reversing the direction of surface activation and 
swimming speed). In such a case we must have that the pitching angle reaches a value — Q 
when the swimmer has returned to the distance h = 2 on its journey away from the sur- 
face. Having determined in this case that a critical initial pitching angle for wall-impact 
is given by 0q ~ —0.4, it might then have been predicted correctly that all wall-impacting 
trajectories in this case lead to final pitching angles of Of ~ 0.4. The near- wall contact 
behavior simply acts to remove any information about the initial pitching angle until the 
trajectory matches that of the critical case (outside of a horizontal shift along the wall). 

For the spherical squirmer that swims with speed unity in a quiescent fluid we have 
(3 = 1/2 (see 4.18), and the dynamics predicted by the reduced model are set by 

h = i 1 + p) sin(0) ' §= w i 1 + t) cos{e) - (5 - 3) 

Can we predict the scattering behavior described above analytically? In order to derive 
a simple estimate of the swimming behavior, let us linearize the motion about = and 
assume (30 <C h 3 . Then the translational swimming speed is not varied by the wall beyond 
its effect on the swimming angle 0, and we have 

h = 0, 0=^, (5.4) 

with h(0) = ho and 6(0) = h(0) = , where the geometrical dependence has been 
absorbed into the singularity strength, (3 = (l + 3/ 1 /2)/3. Integrating this system, we 
find (for trajectories that do not impact the wall) that 
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or equivalently, 

* ^ (5.6) 



8 \ hl h 3 



Taking t — > oo, we know from the simulations that all squirmer trajectories have h — > Of 
and h — > oo; upon insertion, we find a final pitch angle for non- impacting squirmers of 




9 t = \^ + 7iS- ( 5 - 7 ) 



The distance to the wall h = h w {9) = + e 2 ) — (1 — e 2 ) cos(20)/\/2 here corre- 
sponds to wall impact. First we ask: which initial distances and orientations lead to wall 
impact? Writing the pitching angle at the time of impact as W , the centroid will be 
located at a distance h w (9 w ) at that time. Inserting 9 W and h w into (5.5) gives 



/iq h w (9 v 



(5.8) 



Given the simple expression for the translational velocity h = 9, and that 9 > at all 
times, a wall- impacting squirmer must have 9 W < 0. A curve in (ho, 9o) space separating 
initial conditions for which a squirmer does or does not impact the wall may then be 
deduced by setting 9 W — 0, leaving 

°K ~ P * (59) 

y + 8^-8^' (5 - 9) 

or 




In the case of spherical squirmers, for ho = 2 as in Fig. 9, this evaluates to approximately 
9o = —0.33, which is slightly smaller in magnitude than the critical angle 9 ~ —0.4 found 
in the full simulations. Initial pitching angles smaller than 9q = —0.33 are predicted to 
lead to wall collisions. As the initial position becomes more distant from the wall, the 
critical angle increases in magnitude, and depending on the value of /3 there may not 
exist an initial orientation such that the swimmer impacts the wall. (The special case 
of swimming directly towards the wall and guaranteeing impact, 9 = — n/2, is not 
accounted for in the linearized system.) 

After the wall contact there is rotation for a time T(9 W ) until the pitching angle reaches 
9 = 0, at which point the body is predicted in this approximation to separate from the 
wall and depart to a final pitching angle. The final pitching angle may be determined by 
setting 9o = and ho = e in (5.7), giving 




f = \£*> ( 5 - n ) 



which in the spherical squirmer example returns the value 9f ss 0.35. The predicted 
scattering angles from (5.7) (and 5.11 for 9 n < —0.33) are shown in Fig. 9a as circles. 
The simple estimates derived above provide a remarkably accurate depiction of the full 
interaction dynamics with the surface. 

The role of the body geometry is noted here simply by inserting the exact source dipolc 
strength as a function of aspect-ratio, which is a monotonically increasing function of e 
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~ ' t/2 

Figure 10. Contours of the final pitching angle Of of a squirmer, initially placed at h = 2, 
as a function of the initial pitching angle 6q and the squirmer aspect-ratio e. (a) As predicted 
integrating (5.3) numerically, and (b) as predicted using the linearized approximation described 
in the text. Once again the agreement is reasonable when the initial and final pitching angles are 
not very large. The final pitching angle reached by wall-impacting swimmers is not monotonic 
in the body aspect-ratio. 



from its minimum of for a slender squirmer to 1/2 for a spherical squirmer (see 4.18). 
Scattering behavior depends not only upon the source dipole strength, however, but also 
upon the geometrical factor (1 + 3P/2). In particular, we note that (5.11) may without 
difficulty be written solely as a function of the aspect-ratio, e, and we find that Of is a 
monotonically decreasing function on the domain e € [0, 1]. As e — > 1 (a slender rod), 
we have 9f — > \/2/A. The wall interaction is predicted to be strongest for a spherical 
squirmer, and the final angle at which the body swims away from the wall is predicted 
to be greatest in this case. 

Using the far-field theory we may rapidly produce a contour plot of the final pitching 
angle as a function of both the body aspect-ratio e and initial pitching angle 0q , where 
again we initially set the body at a distance h = 2 from the wall. Figure 10a shows 
the predicted values obtained by integrating (5.3) numerically, while Fig. 10b shows the 
analytically predicted values from the simplified approach above. Once again, as in Fig. 8, 
the agreement is quite good when the initial and final pitching angles are not very large. 
The numerical integration of the complete far-field theory indicates that the final pitching 
angle reached by wall-impacting swimmers is not monotonic in the body aspect-ratio. 
This effect is not captured in the linearized approximations above. Making predictions 
about the consequences of wall-impact behavior for non-spherical bodies is complicated 
by the added repulsion force in the numerical simulations. In particular, the repulsive 
force breaks a time symmetry in the otherwise time-reversible structure of the Stokes 
equations, so arguments depending on this symmetry are generally invalid. Nevertheless, 
when the initial pitching angle is not too large in magnitude the hydrodynamics are often 
sufficient to keep the swimmer sufficiently far from the wall, and in such situations the 
repulsive force never plays a role. In these cases we find excellent agreement between 
predictions and the full simulation results. 

Finally, when the body is in contact with the wall the h and 9 dynamics are effectively 
decoupled. Assuming that during wall contact h = h w (9), and expanding 9 [see (5.3)] for 
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small pitching angles, we may integrate 



3£ 
8? 



1 



e 2 + o(0 4 ), 



to find 



0(f) = — tanh 

9e 



8e 4 



+ tanh 1 (g e 0„,) 



(5.12) 



(5.13) 



where q e = y/2 — 3e 2 /2/e, and the body is taken to impact the wall at t = with 
pitching angle W . The time required for the body to rotate to = (at which point the 
body departs from the wall) is then 



T(6 W ) = tanh" 1 {-q e 6 w ) . 

3<?e/3 



(5.14) 



6. Bacteria-like polar swimmers: geometric asymmetry 

6.1. Model swimmers and computational adjustments 

In order to probe the role of polar body geometry in hydrodynamic interactions with 
walls for more biologically relevant swimmers we will enlist the help of a different class of 
model organisms, as illustrated in Fig. 11. The body is composed of an inert spheroidal 
head with centroid x a (on which the no-slip condition is applied) and an active propelling 
rod of dimensionless length L. The system is made dimensionless by scaling on the semi- 
major axis length of the head and on the free-space swimming speed as before. The 
propulsive mechanism is a prescribed uniform distribution of a tangential force per unit 
length on the rod, e • f(s) = — J 7 , where s G [0, L] is the arc-length parameter along the 
rod and T is determined numerically so that the body swims with unit speed in free- 
space. The fluid velocity on the rod is thus composed of a rigid body motion in addition to 
a tangential slip velocity along the long axis, u s (s)e, which must be determined through 
force and torque balance. 

While a uniform distribution of force along the rod is chosen for the sake of simplicity, 
variations in force distribution are typical in planar flagellar undulations and can play an 
important role in setting swimming trajectories. For example, Smith & Blake (2009) have 
shown that the wavenumber in spermatozoan swimming influences the force distribution 
along the flagcllum, which in turn affects the pitching dynamics of the organism near a 
surface. Variations in the force distribution can in fact be an order of magnitude larger 
than their mean values (see Johnson & Brokaw 1979; Smith & Blake 2009). To what 
extent the change in pitching dynamics with wavenumber is due to a change in the time- 
averaged singularity strengths, or due to the time-variation of those strengths, remains 
to be seen. It is possible that the assumption of uniform force distribution is more suited 
to studying bacteria with helical flagella, for instance (see Lighthill 1996). 

Numerically, the rod is represented simply as a line of Stokeslets, which may be included 
into (4.8) without difficulty. The spacing between the points on the rod is chosen to be 
the same mean distance as between points on the body, h gr id- In order to compute the 
singularity strengths, the integrals (4.14)-(4.16) are still appropriate for determining the 
contribution of the head to the singularity strengths (and the u integrals disappear due 
to the no-slip condition that is applied on the spheroid). But to these values we add 
contributions from the rod as follows. Again aligning the body along e = x for the sake 
of presentation, the singularities emanating from the rod center have strengths (from 
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Figure 11. A polar bacteria-like model swimmer directed along e composed of a spheroidal 
head with aspect-ratio e and an actively pushing rod. (The system is still made dimensionless 
by scaling lengths upon the semi-major axis of the head). 



(4.14)-(4.16)): 



XT J 8=0 



P' = 0, 



7=7^/ (x(s)-X) 2 f x (s)ds, 
16n J 8=0 



(6.1) 

(6.2) 
(6.3) 
(6.4) 



where A = —(1 + L/2), x(s) = —(1 + s), f x = x ■ f, and we have included F' as the 
coefficient of the Stokeslet singularity (which balances the Stokeslet singularity strength 
of the moving sphere due to the zero net force condition) . In order to express the singu- 
larities as emanating from the spheroid center, x , the singularities above may simply 
be shifted, yielding the natural expressions 



a 



XF' 



1 



x(s)f x (s), 



Prod — 0, 



Irod = l' + Aa' + 



A 



1 



F' = — / xisffJs). 
2 16tt J s=0 W J K ' 



(6.5) 
(6.6) 
(6.7) 



Adding these last integrals to the expressions in (4. 14)- (4. 16) returns the singularity 
strengths of the entire spheroid-rod swimmer. Similar expressions for the singularity 
strengths generated by a body composed of multiple ellipsoids may be produced in a sim- 
ilar way, shifting the dipole and quadrupole strengths as required, suggesting a tractable 
method for constructing the singularity strengths for an arbitrary swimmer geometry. 
The slenderness of the rod allows us to neglect the slip velocity in the determination 
of the singularity strengths here. The singularity strengths so computed using the full 
boundary integral solution to the Stokes equations (4.8) are shown in Fig. 12 for a range of 
rod lengths L and body aspect-ratios e. One simple estimate of the singularity strengths 
for any body aspect-ratio may be deduced by assuming no hydrodynamic interactions 
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a (Stokeslet dipole) 



(Source dipole) 



7 (Stokeslet quadmpole) 




Figure 12. Computed singularity strengths (about the sphere centroid a; ) for sphere-rod swim- 
mers of varying rod length L and head aspect-ratio e. The quadrupole strength is always nega- 
tive, due only to the placement of the singularities at x a . 



between the spheroid and rod, yielding 



1 
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y 



7 



1 



(0.8) 



The Stokeslet quadrupole strength is always negative, which seems to be at odds with 
our exploration of the ellipsoidal Janus swimmers. The discrepancy is due to our choice 
to write the singularities as emanating from the center of the spheroid, x . By shifting 
the singularities to another point along the body, such as the center of reaction, the 
Stokeslet quadrupole strength can take on positive or negative values. 

0.2. A modified Faxen Law for spheroid-rod swimmers 

Faxen's Law as described by (2.18)-(2.19) is inadequate for the body geometry in Fig. 11. 
While it is common practice to represent microorganisms of arbitrary geometry and 
propulsion as prolate ellipsoids for the sake of mathematical tractability, this approach 
in general can lead even to qualitatively incorrect predictions (to be more precise, quali- 
tative departures from the values found in the full simulations). To derive a more appro- 
priate representation of the effect of an ambient flow on the spheroid-rod model swimmer, 
we apply the usual Faxen Law to the sphere, but couple this to a slender body represen- 
tation of the rod to describe its interaction with the surrounding fluid. For simplicity we 
neglect hydrodynamic interactions between the two, but require the spheroid-rod system 
to move as a rigid body while satisfying the zero net force and zero net torque conditions. 
This approach returns surprisingly good agreement with a "perfect" Faxen Law for our 
purposes, in which the full boundary integral formulation is used to solve for the rigid 
body motion in a prescribed ambient flow. The errors in what follows are due instead 
primarily to the breakdown of the far-field representation of the wall-induced flow field. 

The spheroid is taken to move through the fluid with dimensionless translational ve- 
locity U and rotational velocity fl. In this approximation, the motion of the spheroid 
then applies the following (dimensionless) force and torque to the surrounding fluid: 



and 



F s = Ott [X A ee' + Y A {I - ee')] {U - u*{x a )) , (0.9) 



T s = 8tt [X c ee' + Y c (I - ee')] [Q - u*) - Si:Y H e x E*(x a ) ■ e, (0.10) 



where X A = (8/3) X 3 /(-2 X + (1 + X 2 )L X ), Y A = (10/3) X 3 /(2x + (3% 2 - 1)L X ), X c 
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(4/3) X 3 (l - X 2 )/(2 X - (1 - X 2 )L X ), Y c = (4/3) X 3 (2 - X 2 )/(-2x + (1 + X 2 )L X ), Y H = 
(4/3) X 5 /(-2x + (l + X 2 )i x ), L x = log[(l + x)/(l- X )], and X = (Kim & Karrila 



1991). Recall that u* in our case is the fluid velocity created by the image singularities 
as described in §11. 

Meanwhile, a local slender body approximation is used to describe the fluid force 
imposed by the motion of the rod, in which the force per unit length is written as 



where c = log(4/e 2 e), e is the aspect-ratio of the slender rod, and we estimate the 
slenderness of the rod by taking e = Sy/e/(2L) (see Gray & Hancock 1955; Childress 
1981). Since the rod moves with rigid body motion, u(s) = U — s Q x e. Note that 
the slip velocity does not enter into the Faxen-type approach here. Also, we refer the 
reader to Bouzarth & Minion (2011), who explore slender body theory using regularized 
Stokcslets in much greater detail. 

The motion of the rod applies the following force and torque (about x Q ) on the sur- 
rounding fluid: 



The zero net force and zero net torque conditions, Fg + Fr = 0, Tg + Tr = 0, produce 
a linear system which may be solved for the resulting translational velocity U and rota- 
tional velocity fl. Into these expressions, we simply insert the singularity image systems 
for the tilted force dipole, source dipole, and force quadrupole as before, as it*. 



We are now equipped with the tools to compare the full simulations with the far-field 
predictions derived above. Fixing the pitch angle to 9 = — 7r/6, and the head aspect-ratio 
to e = 1, we show in Fig. 13 (top row) that the simple Faxen Law developed above 
captures very well the adjustment to the swimming dynamics when the body is far from 
the wall, even without the inclusion of hydrodynamic interactions between the sphere 
and the rod. Scanning to even smaller distances (Fig. 13, bottom row), we find very 
accurate quantitative similarity between the predictions and the simulation results down 
to approximately h = (L + 2) (the total length of the sphere-rod body). When the body 
is yet closer to the wall, the far-field prediction retains qualitative accuracy for a range 
of distances, but eventually becomes qualitatively incorrect. As seen in the comparison 
of U x and U z on the bottom row of Fig. 13, the translational velocities predicted by 
the far-field theory in this particular case take the wrong sign and yield non-physical 
behavior. For all but these very small distances from the wall, the far-field prediction 
with the Faxen Law of the previous section provides a quick and accurate approach for 
predicting the wall-induced adjustments to the swimmer trajectories. 

Again we have verified by comparison with a "perfect" Faxen Law (using the full 
simulations to determine the effect of the image singularities) that this degradation is 
not due to the breakdown of the modified Faxen Law but due to the breakdown in the 
far-field approximation of the wall effect. 




(6.11) 




(6.12) 



6.3. How accurate is the far-field representation? 
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Figure 13. Comparison of full simulation and far-field predictions for 9 = — tt/6, with computed 
values indicated by symbols and far-field predictions indicated by lines. The top row shows results 
for bodies well separated from the wall, (h — 1)/(L + 2) £ [2,4], while the bottom row shows 
results for bodies very close to the wall, (h— l)/(L + 2) € [0, 2]. The simple Faxen Law developed 
in the text captures well the adjustment to the swimming dynamics when the body is far from 
the wall, even without the inclusion of hydrodynamic interactions between the sphere and the 
rod. For bodies very close to the wall the far-field prediction retains qualitative accuracy for a 
range of distances, but eventually becomes qualitatively incorrect. 

7. Employment of the reduced model: spheroid-rod swimmers 

7.1. Wall-induced rotations 

As in §IV, we now consider the initial effect of the wall on bacteria-like swimmers of 
varying geometry and propulsive activity. Our previous analysis for ellipsoidal model 
swimmers suggests that the ratio of head-size to flagellum length might play a critical 
role in setting the wall-induced rotational behavior. 

We begin by fixing a swimmer with a spherical head to a distance h = 5 from the wall, 
and the pitching angle to 9 = 0, and plot the wall-induced rotation as a function of rod 
length in Fig. 14a. And indeed, we see that there is a critical rod length at which the 
rotation generated by the proximity to the wall reverses sign. The results suggest that for 
organisms with small flagella relative to cell-body size, the rotation will act to redirect the 
swimming away from the wall; similarly, for organisms with longer flagella relative to the 
cell-body size, the wall-effect will be to redirect the swimming towards the wall, increasing 
the likelihood of wall impact. This critical rod length was determined by simulations to 
take the value of approximately L* = 4.5, or slightly more than twice the cell-body 
diameter. The far-field prediction also shows this transition in the direction of rotation. 
The prediction is in close quantitative agreement with the full simulations for small rod 
lengths (where the far-field approximation is expected to be more accurate), but slightly 
underestimates the critical rod length for direction reversal, L* a = 3.0. Nevertheless, the 
important qualitative features of the wall-effect are captured by this simple means of 
estimation. 

Now let us place the bacteria-like swimmer much closer to the wall. Fixing h = 1.5, 
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Figure 14. (a) Rotation rate for swimmers with spherical heads and varying rod lengths L 
with 8 = and h — 5 fixed, (b) The same, for a fixed distance h = 1.5. (c) The equilibrium 
pitching angle determined using fully time-integrated simulations (h and 9 are determined by the 
dynamics) as a function of L is shown as a dark line. The pitching angle in the steady swimming 
state is positive for rod lengths L < L* w 2, and negative for rod lengths L > L* . The negative 
equilibrium pitching angle does not vary significantly as the rod length is increased, but takes 
on a minimum value at L as 4. The far-field prediction of equilibrium pitching angle for three 
fixed heights are included for comparison, (d) The critical rod length for which there is rotation 
direction reversal as a function of the distance to the wall, h, as determined using full numerical 
simulations. 



we plot the rotational velocity as determined by simulation and far-field prediction in 
Fig. 14b, again as a function of rod length. While the quantitative error in the far-field 
prediction increases, the simple estimation is still able to capture a qualitative difference 
in pitching dynamics for rod lengths below and above a critical value. The result that 
there exists a critical rod length for rotational direction reversal appears to be robust, 
given that the adjustments to this critical rod length from h = 5 to h = 1.5 are quite 
small. 

Loosely applying the theoretical predictions for ellipsoidal body geometries described 
in §111 (for instance the linearized expression shown in 5.1), we might expect that a 
large source dipole strength relative to the Stokeslet quadrupole strength would lead 
to an induced rotation away from the wall. In drawing this comparison recall that the 
singularity strengths represented in Fig. 12 are centered about the centroid of the head, 
so that the source dipole should be compared instead to an appropriate combination of 
the Stokeslet quadrupole and a shifted Stokeslet dipole. This rough comparison provides 
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intuition as to why an increase in rod length for a given head shape, which increases the 
Stokeslet dipole and quadrupole strengths more rapidly than the source dipole strength, 
can lead to a reversal in the wall-induced rotation. 

In general the rotational effects described above are coupled to the wall-induced at- 
traction or repulsion in setting the complete swimming trajectory of a given organism. 
We already see, however, that the geometry of a swimming microorganism is very likely 
to play an important role in such dynamics as wall-entrapment or adhesion and the de- 
velopment of biofilms. While organisms with small flagella relative to cell-body size may 
still be hydrodynamically (or otherwise) pinned to a solid surface, their orientation in 
such a state can be directed away from the wall and the entrapment may be unstable to 
sufficient perturbation. Meanwhile, for microorganisms with longer flagella, the rotational 
component of the wall-effect is likely to act in tandem with hydrodynamic attraction to 
greatly increase the likelihood of wall-impact and long-term entrapment. 

7.2. Equilibrium pitching angles 

We now allow the polar bacteria-like body to swim until it finds an equilibrium pitching 
angle. From the previous consideration we expect to find a critical rod length below which 
the equilibrium angle is positive, and above which the equilibrium angle is negative. The 
bacteria-like swimmers studied here have a > 0, so that the effect of the wall at leading 
order (when 9 <~ 0) is to attract the swimmer to the surface. For a small enough initial 
distance h(t = 0) we find that these swimmers (with L <G [0.5,7]) all become pinned to 
the wall. For larger initial distances, not all swimmers will impact the wall and find such 
an equilibrium state. Particularly, for organisms with small flagella relative to the cell 
body size, there will be a competition between the rotation away from the wall (which 
redirects the swimming velocity) and the hydrodynamic attraction towards the wall. 

Figure 14c shows the equilibrium pitching angle determined using fully time-integrated 
simulations (in which ft and 9 evolve dynamically) as a solid line. As in §V, a repulsion 
force is included in order to keep the body from making contact with the wall of the form 
F rep = Ae~ Bd /{l — e~ Bd ]z, where d is the minimum distance between the body surface 
and the wall, and here we take A = 100, B = 100; with this repulsion the body does not 
come closer than ft as 1.01 for any rod length. Starting from various initial distances (close 
to the wall) and pitching orientations, and evolving in time until a steady equilibrium 
is reached, the plotted equilibria appear to be the only stable equilibria. The pitching 
angle in the steady swimming state is seen to be positive for rod lengths L < L* ss 2, 
and negative for rod lengths L > L* . The negative equilibrium pitching angle does not 
vary significantly as the rod length is increased, but takes on a minimum value at L ss 4. 

We include in Fig. 14c the equilibrium angles predicted by the far-field theory for 
fixed distances ft = 5, ft = 1.5, and ft = 1.1. The far-field theory provides a reasonable 
qualitative prediction of the trapped equilibrium state, but only when the distance to the 
wall is not taken to be too small. The far-field theory breaks down when the swimmer 
is close to the wall, and the pitched-down equilibrium states are not predicted in this 
range of rod lengths. Higher order terms in the multipole expansion not included here 
(or lubrication effects) are thus found to play a vital role in determining the very-close- 
contact behavior, just as was generally the case for the model Janus swimmers. 

Note that the critical rod length for direction reversal does depend on the distance ft. 
At a distance ft = 5 (fixing = 0), bodies with rods longer than L* = 4.5 will rotate 
downward towards the wall, while only bodies with rod length greater than L* = 5.1 
rotate towards the wall when ft = 1.5 (as determined using the full simulations; see 
Figs. 14a-b). When the body is in very close contact with the wall, ft = 1.1, then the 
critical rod length is approximately L* = 2.0 (see Fig. 14c). The non-monotonicity of the 
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critical rod length for rotation direction reversal is shown in Fig. 14d (on a semi-log scale 
in the inset), and this critical length is found to decrease to zero as the spherical head is 
moved closer to the wall. The motion of a sphere moving close to a plane wall has received 
classical attention (see for instance Goldman, Cox, & Brenner 1966), where it has been 
shown that the effect of lubrication forces act to "roll" a sedimenting sphere along the 
surface. Correspondingly in the case studied here, the lubrication forces are expected to 
contribute to a downward pitching motion, and thus to reduce the critical rod length 
for direction reversal at very small distances. This is not to suggest that, absent a steric 
repulsion, all swimmers would eventually turn towards the wall to a negative equilibrium 
pitching angle. A swimming velocity directed away from the wall can be sufficient to 
balance hydrodynamic attraction at a distance h such that lubrication forces are not 
dominant. 

Finally, another geometrical factor of possible importance in wall-induced swimming 
behavior is the shape of the cell-body. For our last consideration, let us return to a 
regime where the far-field theory retains a modest predictive power. Figure 15a shows 
the equilibrium pitching angle for a range of rod lengths L and head aspect-ratios e 
when the head is fixed at a distance h = 5 from the wall, as determined through full 
numerical simulation. The direction of rotation for bodies swimming parallel to the wall 
at this distance may be inferred by continuity. For a given head shape, we find that there 
is a critical rod length in the range L e (3,5) for which the equilibrium pitching angle 
changes sign, denoted by a dashed line. The equilibrium angles predicted by the far- 
field theory are shown in Fig. 15b. The far-field theory captures the change in direction 
reversal, though underestimates the critical rod length, and also generally overestimates 
the magnitude of the equilibrium pitching angle. The far-field theory does capture the 
decrease in the critical rod length for head shapes of smaller aspect-ratio. 

Of particular note, Figs. 15a-b suggest that organisms with sufficiently long flagella may 
be passively rotated towards the wall regardless of the shape of the cell-body. Similarly, if 
the flagella is sufficiently short the induced rotation may be oriented away from the wall 
for any head shape. Once again applying loosely the exact theory for ellipsoidal geometries 
from §111, as the cell-body becomes more prolate in shape (and hence decreases in volume 
given our choice of nondimensionalization), the source dipole strength rapidly decreases 
in magnitude to zero. This would suggest that the contribution to rotation that is oriented 
away from the wall diminishes with decreasing aspect-ratio, so that the critical rod length 
for direction reversal as e — > should also decrease, consistent with our findings. 

8. Discussion 

In this paper, we have explored a multipole representation of swimming organisms 
to improve our understanding of swimming behaviors near a surface from a general 
hydrodynamic perspective. The simplified model has shown itself to be very useful in 
understanding the behavior of a selection of model swimmers of varying geometry and 
propulsive activity. Our results present some hope that by determining only a small 
number of parameters describing a given organism, much of its hydrodynamic interaction 
with a solid or free surface, or with another similar organism, might be predicted without 
the need for large scale numerical simulations. The framework described here may be of 
use in understanding as yet puzzling behaviors in near-wall microorganism behavior, 
and conversely may be of great use in the design of moving structures to perform tasks 
near a surface. One of the most general findings in our study has been that swimming 
bodies with small inert heads and long active propulsive mechanisms are those bodies 
which rotate nose down towards the wall due to hydrodynamic interactions, while other 
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Figure 15. (a) Contours of equilibrium pitching angle 8 for fixed height h = 5 computed using 
the full simulations, (b) The same, but as predicted using the far-field approximation. As in 
the case of the Janus swimmers, bodies with small no-slip regions (cell bodies) and long active 
regions (flagella) are pitched downwards passively by the hydrodynamic interaction with the 
wall, and vice versa for large cell bodies with small flagella. 



swimmers rotate nose up away from the wall. Among other reasons, the question of 
equilibrium geometry (i.e. orientation) is important because it impacts the response of 
a wall-bound organism to a shear flow, a topic of intended future study along the lines 
considered herein (Hill et al. 2007; Kaya & Koser 2009). 

Avenues of interest for future study include time-dependent singularity strengths (due 
perhaps to periodic undulations in propulsion mechanism) and their consequences on 
wall-induced behaviors and trajectories. Such variations may well change the nature of 
the equilibrium states near the wall, and this is suggested by the work of Shum, Gaffney & 
Smith (2010) and Elgeti, Kaupp & Gompper (2010). The matching of far-field predictions 
with near-field lubrication effects may also prove fruitful in certain cases, which would 
allow for a more complete swimming trajectory and hydrodynamic entrapment analysis. 
The framework described above might be pursued further to predict long time trajectories 
for a number of swimmers. 

While not considered here, we have included the stress-free boundary condition veloci- 
ties in Appendix B, and there arc likely many interesting phenomena which can be tuned 
or explained by varying the boundary conditions at the surface in precisely this way. The 
qualitative behaviors predicted by the far-field consideration for the Stokeslet dipole, 
source dipole, and Stokeslet quadrupole do not vary dramatically for such boundary con- 
ditions, but the same cannot be stated for the rotlet dipole singularity (as described in 

§ni). 

Among other physical and biological effects, we have ignored here the consequences of 
Brownian motion in order to take a clearer look at hydrodynamic effects. The interac- 
tion of both hydrodynamic and Brownian effects in swimming microorganisms has seen 
significant recent interest. Li, Tarn & Tang (2008) and Li & Tang (2009) have shown 
that wall collisions and rotational Brownian motion alone can give rise to accumula- 
tion of microswimmers near a surface. Reintroducing the hydrodynamic effects of the 
wall, the same authors also show the ways in which Brownian motion is amplified by 
hydrodynamic interactions with the wall. Armitage & Macnab (1987) have shown that 
Rhodobacter sphaeroides is apparently governed significantly by rotational Brownian mo- 
tion alone. Equipped with a more complete theory for the wall effects on motility, we 



Hydrodynamics of self-propulsion near a boundary 



33 



expect that the interplay between hydrodynamics and stochastic processes may be easier 
to distinguish in many cases. While Brownian motion and steric effects have been argued 
to drown out fluid mechanical interactions in swimming bacteria (and bacteria-wall in- 
teractions) by Drescher et al. (2011), this is less likely to be the case for organisms of 
larger scale. 

We acknowledge the support of the NSF through Grant No. CBET-0746285. 



Appendix A. Vector notation for the fundamental singularities 

The fundamental singularities of the Stokes equations have a long history in the study 
of fluid dynamics, and their tensorial representations may be found in Pozrikidis (1992). 
We refer the reader also to the work of Chwang & Wu (1975). Here we include the 
vector notation of the fundamental singularities for reference. The Stokeslet, and Stokeslet 
dipole, quadrupolc, and octupole may be written as 




(Al) 



_ . , . 1 / (d ■ x)e — (e • x)d — (d ■ e)x 3(e • x)(d ■ x)x\ 
G D (x- d, e = 1 — 12 ( ^ '- ^—-^ * '— + -i '— , 

\ x \ V \ x \ \ x \ J 




(A 2) 



G Q (x; c, d, e) = — {{d ■ e)c + (c • e)d - (c • d)e - + 



)■ 



(A3) 



G (x;b,c,d,e) = r—^ 




(A 4) 



where 



.91 = [(d ■ e)(c • x) + (c • e)(d ■ x) + (c • d)(e ■ x)] x 

+ (d ■ x)(e ■ x)c + (c • x)(e ■ x)d — (c • x)(d ■ x)e, 



(A 5) 



.92 =[( c • d)(e ■ x) + (c • e)(d ■ x) + (d ■ e)(c • x)]b 
+ [(d • e)(b ■ x) + (b ■ e)(d • x) + (b ■ d)(e • x)]c 
+ [(c ■ e)(6 • x) + (b ■ e)(c ■ x) + (b ■ c)(e • x)]d 
- [(c • d)(b -x) + (b- d)(c -x) + (b- c)(d ■ x)]e 
+ [(c • d)(b • e) + (c • e)(6 ■ d) + (b ■ c)(d ■ e)]x, 



(A 6) 



#3 =(c ■ x)(d ■ x)(e ■ x)b + (b ■ x)(d ■ x)(e ■ x)c 
+ (b ■ x)(c ■ x)(e ■ x)d — (b ■ x)(c ■ x)(d ■ x)e 

+ (b ■ e)(c • x)(d ■ x) + (6 • d)(c • x)(e ■ x) + (b ■ c)(d • x)(e ■ x) 



+ (c ■ d)(b ■ x)(e ■ x) + (c • e)(6 ■ x)(d ■ x) + (d • e)(6 • aj)(c • x) 
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Figure 16. A model swimmer near a wall. The body swims in a direction e, which may be 
decomposed into a horizontal part along and a vertical part along z. The pitching angle 
9 = cos _1 (e • e^,) is indicated, as is the distance from the wall h of the body centroid. cj> denotes 
the angle of the horizontal swimming component relative to the x axis, and for the present study 
any variations in cj> are decoupled from the pitching and height dynamics of 9 and h. Ignoring 
the effect of the rotlet dipole, we may set (f> — and = x without loss of generality. 



Meanwhile, the source doublet and higher order source singularities in vector form are 
s 1 ( 3 (as • e)x\ , , . 

Q{x , d, e) = f£ ( (d.x)e + (e.x)d+(d.e)x 5(e ■ x)(d ■ x)x \ 

3 / 5<74 35(c • x)(d ■ x)(e ■ x)x 



0(a;;c ) d,e) = r -| ¥ f(c-d)e+(c-e)d+(d-w^ , 



(A 10) 



where 

54 = (d • x)(e ■ x)c + (c ■ x)(e ■ x)d + (c ■ x)(d ■ x)e 
+ [(d • e)(c • x)x + (c • d)(e ■ x) + (c • e)(d • a;)] x. (A 11) 

Finally, the rotlet and rotlet dipole may be written as 

R(x;e) = ^^, (A 12) 

„ , , . dxe 3(d-x)(exx) ,,..„.. 
i? D (:E; d, e) = + ^ ^ ^- (A 13) 



Appendix B. "Tilted" image systems of the fundamental singularities 

The image systems for the fundamental singularities (those axisymmetric about the 
swimming direction e) are dependent upon the distance of the singularity to the wall, 
h, and the orientation, e. The orientation e is written in polar and azimuthal angles 
describing the swimming orientation via e = cos(6>) cos(4>)x + cos(6*) sin(0)y + sin(0)z, as 
shown in Fig. 16. We proceed to show the effect of a no-slip wall and a stress-free surface 
on a body of aspect-ratio e for each of the fundamental singularities considered in (2.17). 
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B.l. Tilted Stokeslet image system 

The image system for a Stokeslet near a wall is well known (Blake & Chwang 1974). The 
image system for a Stokeslet directed along e = cos(#) cos(<p)x + cos(#) sin(</>)y + sin(#)£ 
in the fluid is simply a linear combination of these known image systems, and we write 

G*(x - x* a , e) =cos(0)( - G(e ) + 2hG D ( e4> ,z) - 2h 2 D(e^ 

+ sin(0) ( - G(z) - 2hG D {z, z) + 2h 2 D(z)^j , (B 1) 

where = cos(<f))x + sin(0)y. For the image systems described here and below, we use 
the shorthand notation G(e) = G(x — x* a ;e) (the singularities are placed at the image 
point inside the wall). 

The addition of the Stokeslet singularity and the image set above placed at the image 
location x* renders the fluid velocity zero on the wall exactly for arbitrary distances h 
from the wall, 



u(a:)U=o= J- (G(x-x ;e) + G* (x - x* ;ej) = 0, (B2) 

OTT/I V / z=0 

where x a * = x a — 2hz is the position of the image singularity (inside the wall). Faxen's 
Law (2.18)-(2.19) is now used to determine the wall-induced translational (u) and rota- 
tional (17) velocities of a (prolate ellipsoidal) body with a Stokeslet placed at its centroid 
z ■ x a = h: 

3cos(g) 3sin(fl) - 3rcos(^(l + sin 2 (g)) 
= jj— , u z = V= — 2 (B3) 

Here we have written u = u^e^ + u z z and 17 = (z x e^) (so that 9 — —Qj,± ). There 
is no rotation outside of the z x component, which is assured by symmetry. 

It is simpler to determine the velocities induced by the presence of another swimmer 
symmetrically placed opposite an imaginary plane, or equivalently by the presence of a 
stress-free boundary, at z = 0; we need only apply an image system which consists solely 
of the mirror image of the primary singularity, or G(x — x* a ; e*), where e* = e— 2(e • z)z. 
These induced velocities are 

. cos(0) . sin(0) - cos(g)(l + 3rsin 2 (g)) 

M ^^T' Uz = — fT' Q ^ = 47,2 ( ) 

A Stokeslet aligned parallel to a stress-free boundary therefore increases the wall-parallel 
velocity component, while the same Stokeslet induces a decrease in the wall-parallel 
velocity near a no-slip wall. 

B.2. Tilted Stokeslet dipole image system 

The image system for a tilted symmetric Stokeslet dipole near a wall is not simply a linear 
combination of wall-parallel and wall-perpendicular Stokeslet dipole image systems, as 
noted by Blake & Chwang (1974). Instead, we may determine the image system for a tilted 
Stokeslet dipole by taking derivatives of the Stokeslet image system, where derivatives 
along z incorporate the h-dependence of the singularity coefficients. The resulting image 
system is 

G* D (x - x* a ; e, e) = cos 2 (6>) (-G D (e f e ) + 2hG Q (e <tn e , z) - 2h 2 Q(e (j} , e )) 

+ sm 2 (6) (-G D (z, z) + AhD{z) + 2hG Q (z, z, z) - 2h 2 Q(z, z)) 

+ sin(0) cos(0) (G D (e , z) + G D (z, e ) - 4hD( e<t> ) - 4/iS Q (e^, z, z) + Ah 2 Q(e^ z)) 

(B5) 
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Faxen's Law returns the wall-induced translational and rotational velocities now of 



3sin(26») „ 



3(1 - 3sin 2 (<9)) 



3sin(20) 
16h 3 



l + ^(l + sin 2 (0))). 



(B6) 



Once again, as for all the symmetrized singularities about e, there is no rotation outside 
of the z x e0 component, and we have written f2 = Q^x (z x e^). For a swimming body 
with a positive Stokeslet dipole component (a pusher), for small angles 9, the wall induces 
an attraction and a rotation acting to align the body with its long axis parallel to the 
wall. 

For induced velocities near a mirror-image swimmer or near a stress-free surface we 
instead find the qualitatively similar effects: 

fi* = 0, ^ = -^(l-3sin 2 (0)), ^ = ^M(l + rsin 2 (0)). (B 7) 

B.3. Tilted Stokeslet quadrupole image system 

Similarly, directional derivatives of the tilted Stokeslet dipole system yield the image of 
the axisymmetric Stokeslet quadrupole: 

Gq(x - x* a ; e, e) = cos 3 (6») (-Gq^, e^, e^) + 2hG (e < j > , e^, z) - 2h 2 0(e 4n e^, e^)) 
+ sin 3 (0) (3G Q (z, z, z) + 4D(z) - 2hG (z, z, z, z) - 8hQ(z, z) + 2h 2 0(z, z, zj) 

+ cos 2 (#) sin(0) 3G Q (e0, e^, z) + 2G Q {z, e , e^) - 8hQ(e tj> , e^) 
- 6/iG o (e0, e^, z, z) + &h 2 0{e^, e<p, z) 

+ sin 2 (fl) cos(0) - Gq(z, z, e ) - 6G Q (e , z, z) - 4D(e ) + 6/iG o (e , z, z, z) 



16/iQ(e ,z)-6/i 2 O(e ,z,z) 



(Bi 



The wall-induced translational and rotational velocities associated with the tilted 
Stokeslet quadrupole are surprisingly simple: 

cos(#) , „_ . o.-.v . sin(#) 



4/j 3 



(7-9sin 2 (6»)) 



3cos(6>) 
8h 4 



1 - 3sin 2 (0) + j (11 - 3sin 4 (fl))^ . 



(B9) 
(BIO) 



The induced velocities near a mirror-image swimmer or near a stress-free surface are 



cos(0) , . 2//) x\ 



sin(6>) 



sin 2 (0) 



V = (! - 5sin 2 (0) + r (2 + sin 2 (0) - 3sin»)) . 



(Bll) 
(B12) 



B.4. Tilted source dipole image system 

As in the case of the tilted Stokeslet singularity, the tilted dipole image system is simply 
a linear combination of wall-parallel and wall-perpendicular image systems (see Blake & 
Chwang (1974)), 



D* (x - x* a ; e) = cos(9) (.D(e^) + 2G Q (e <j> , z, z) - 2hQ{e^, z)) 

+ sin(0) ( - 3£>(z) - 2G Q (z, z, z) + 2hQ(z, z)) . 



(B13) 
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The wall-induced translational and rotational velocities are 

cos(6») . sin(0) ~ 3cos(6>) / 3.T . . 2 ._\ . „, 

The induced velocities near a mirror-image swimmer or near a stress-free surface are 
cos(6») . sin(0) ~ 3rcos(9). . 2 ,„.. 

B.5. Tilted rotlet image system 
The image system for a tilted rotlet is given by 

R*(x - x* a ; e) = cos(d) {-R(e^) + 2hD(ej) - G D (ej,z) - G D (z, ej)) - sm(6)R(z), 

(Bl6) 

where = — sin((fi)x + cos((f>)y. The wall-induced velocities on the body are 

u = 0, Q = (5 - 3rsin 2 (0)) * ^ (2 + 3rcos 2 (0)) z. (B 17) 

The image singularity for motion near a mirror-image swimmer or near a stress-free 
surface are qualitatively different than those required near the no-slip wall. The wall- 
parallel component of the image system is simply 

{I - zz')R*{x-x* ,e) = {I - zz')R{-e + 2(e-z)z). (B 18) 

Meanwhile, the mirror-image of the singularity that is perpendicular to the boundary is 
exactly the same singularity, but shifted to the image point, 

z ■ R*(x - x* a ;e) = z ■ R(e). (B 19) 

Combining the above, the induced velocities are in this case 

u=^*\ ft = ^( 1 + 3rsin^))« + |^(2-3rcos»)i. (B20) 

Analagous to the qualitative change for a wall-parallel Stokeslet near a free-slip or 
no-slip surface, here there is an increase in the rotation rate for a rotating body near a 
stress-free surface, and a decrease in the rotation rate for a rotating body near a no-slip 
wall. Interestingly, there are angles for which the direction of rotational precession 
about the z axis depends upon the body aspect-ratio. 

B.6. Tilted rotlet doublet image system 
The image system for a tilted (axisymmetric) rotlet is given by 

R* D (x - x* a ; e, e) = cos 2 (6<) (-R D (e 4> , e ) + 2hQ{e <t> , e£) - G Q (e , ej, z) - G Q (e , z, e 
+ sin 2 (6)R D (z,z), 

+ cos(0) sin(0) R D (z,e^) - R D (e (t> ,z) 

+ 2D(ej) + G Q (e^,z, z) + G Q (ej,ej,z) - 2hQ(ej,z)] . (B 21) 

The induced velocities on the body are 

u = 0, (B22a) 
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" = (6 _ r(1 + 3sin2(e)) * _ si (1 _ 3sin2(0) _ rcos2 w + 3sin2 ^)) £ - 

(B22&) 

The image singularity for motion near a mirror-image swimmer or near a stress-free 
surface is given by 



R* D (x - x* a , e, e) = R D (e - 2(e • z)i, -e + 2(e • z)£). 
The induced velocities are in this case 

3sin(2<9) 1A _ 



u = — - 



16/i 3 



(B23) 



(B24) 



" = - ^^T- ( 2 + r (l + 3sin 2 (0)) + ^4 (1 - 3 sin 2 (0) + rcos 2 (0)(l + 3sin 2 (0))) £. 

(B25) 



B.7. Stresslet image system 

The image system of the Stresslet singularity which appears in the double layer integral 
in (4.2), may be deduced by decomposing it into the fundamental singularities considered 
above. Consider the d, e-Stresslet, as dchncd in (2.15), 



T(x;d, e) 



-6(d • x)(e ■ x)x 



x 



(B26) 



Decomposing the image singularities into operators acting on the components of e, we 
write 



T*(x-d,e) = T* 



e x + T* 



e y +T* 



(B27) 



where 



= -2d x U + Ah d z D(x) + Ah d x D(z) 

+ 2d x G D (x, x) - Ad x G D {z, z) + d y (g d {x, y) + G D (y, x)j 

- d z (g d (x, z) + G D (z, xfj + Ah 2 d x Q(x, x) + Ah 2 d y Q(x, y) 

- Ah 2 d z Q(x,z) + Ahd z G Q (x,z,z) - Ahd x Gq(x,x, z) - Ahd v G Q (x,y,z), 



(B 28) 



= -2d y U + Ah d z D(y) + Ah d y D{z) 

+ 2d y G D {y,y) - Ad y G D {z, z) + d x (g d (x, y) + G D (y, x)j 



d z [G D {y, z) + G D (z, y)J + Ah 2 d y Q(y, y) + Ah 2 d x Q{x, y) 
Ah 2 d z Q(y,z) +Ahd z G Q (y,z,z)- Ahd y G Q (y,y,z) -Ahd x G Q (x,y,z), 



(B29) 
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= -2d z U + 4hd x D(x) + 4hd y D(y) -4hd z D(z) 

- 2d z G D (z, z) - d x (G D (x, z) + G D (z, i)) - d y (G D (y, z) + G D {z, y)) 



+ Ah 2 d z Q(z, z) - Ah 2 d x Q(x, z) - Ah 2 d y Q(y, z) 

+ 4hd x G Q (x, z, z) + 4hd y G Q (y, z, z) - 4hd z G Q (z, z, z). (B 30) 

Recall that all the singularities above are placed at the image point x = x* a , as in Ap- 
pendix B. 

When the Stresslet appears in the boundary integral representation in (4.2), we have 
d = n(y) and e = u(y) — u(x). Note that the single layer integral in (4.2) may be 
eliminated with the introduction of an auxiliary Stresslet singularity density q(x) (Power 
& Miranda 1987; Pozrikidis 1992). In this case, d = n(y) is still determined by the 
geometry but e = q{y) — q(x) is unknown. The framework above may then be used 
to construct the linear operator acting on q, an approach amenable to highly accurate 
numerical schemes which may be made spatially adaptive at no theoretical cost, unlike 
approaches based on the single layer integral which are Fredholm equations of the first 
kind and thus lack the foundational theory underlying the suggested method. 
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